4888 




J. Phys. Chem. 1996, 100, 4888-4892 




Low-Temperature Rate Coefficients of CzH with CH 4 and CD« from 154 to 359 K 


o 


Brian J. Opansky and Stephen R. Leone* + 0 3 & t 

J1IA. National Ins, i, me of Standards and Technology and Uninersdy of, Colorado, Department of Chemistry 
and Biochemistry, University of Colorado, Boulder, Colorado 80309-0440 

Received- Number 6, 1995; In Final Form: January 2, 19915 s 


NASA-CR-204566 


Rate coefficientsfOT^e reactio_n 

^producTby pulsed laser photolysis of ^H^vmable 'rL^effidenteforlhe 

center laser probes the transient removal of C 2 H (X 2 (0,0,0)) m apsorpuon. me , 4 _.cq 

reactions of C,H with ± 1.8) x 

m"^exp[(-650 ± X 61)^]cm’ molecule - ' s"', respectively. The reaction of CtH + CH. exhibits la sigmfican, 
^ ^ ^ t/ A r i, — 2 S 4: 0 2 Temperature dependent rate constants for C2H 4- 

slight negative temperature dependence, which can be expressed as *c : h : - 8.6 x 10 T expl(4/4 ± 
90)/r) cm 3 molecule 's'. 


Introduction 

This work is part of an ongoing project to measure, for the 
first time, low-temperature rate constants of the C 2 H radical 
with various hydrocarbons present in Titan’s atmosphere. In 
planetary atmospheres, such as Titan’s, the ethynyl radical can 
catalyze the dissociation of CH4 to form methyl radicals, CH3. 

C 2 H 2 + hv — C 2 H + H (1) 

C 2 H + CH 4 -C 2 H 2 + CH 3 (2) 

net CH 4 — CH 3 + H ( 3 ) 

The recombination of two methyl radicals produces ethane, 
C 2 H 6 , which can be transported downward to the moon s 
surface. Hence, the reaction C 2 H + CH 4 is of central 
importance in photochemical models of Titan to understand why 
ethane is so abundant in Titan s atmosphere. 

This is the first low-temperature study of the reactions 
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Figure 1. Schematic of the experimental setup. 


these results provide an experimental basis for theoretical studies 
on primary isotope effects. 


C,H + CH 4 (CD 4 ) — C 2 H 2 (C 2 HD) + CH 3 (CD 3 ) (4) 


over the temperature range 154-359 K. A complete study of 
the reaction C 2 H + C 2 H 2 from 170 to 350 K, another critical 
reaction in Titan’s atmosphere, was already reported by this 
laboratory; 3 here, these measurements have also been extended 
t0 J43 — 359 k. Low-temperature rate constants of reaction 4 
will help decide which reaction schemes for ethane production 
are consistent with data from Voyager IRIS (infrared interfer- 
ometer spectrometer) results. 4 - 5 In addition, NASA plans to 
launch a mission called Cassini, which is intended to study 
Saturn and its moon Titan some time in 1996. The rate 
coefficients measured in this experiment will be used to model 
Titan’s atmosphere to compare calculated gas densities with 
measured concentrations from the Cassini mission. In addition, 
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Experimental Technique 

Briefly the kinetics of C2H + CH4 and CD4 are studied using 
transient infrared laser absorption spectroscopy. A schematic 
of the experimental setup is shown in Figure 1. The low- 
temperature kinetics measurements are described in detail in 
an earlier paper. 6 The only significant change was to add 
Brewster windows to the flow cell to incorporate a multipass 
arrangement for the probe laser beam to increase the absorption 

path length. . 

Ethynyl radicals are produced in the meter long variable 
temperature flow cell by a pulsed excimer laser at 193 nm. The 
excimer laser was operated at 55 mJ/pulse at a repetition rate 
of 10 Hz. Acetylene has an absorption cross section of 1.35 x 
10“ 19 cm 2 at 193 nm 7 and a quantum yield of approximately 
0.26 for C 2 H production. 8 After accounting for UV laser loss 
on the windows and absorption in air, the C 2 H concentration 
was calculated to be no greater than 4.6 x 10 cm for the 
highest acetylene pressures (the acetylene number density was 
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Figure 2. (a) Decay of C:H with C 2 H 2 . The solid line is the fitted 
single-exponential decay, (b) Residuals to the fit. The oscillations at 
lower frequencies are due to the noise from the color center laser. 

in the range (1. 6-6.4) x 10 14 cm -3 in the CH_i (CD 4 ) 
experiments). For these experiments a 1000—1500-fold excess 
of CH 4 (CD 4 ) with respect to C2H2 is always present. Contribu- 
tions from secondary or radical— radical reactions can be 
neglected since the time between collisions for C 2 H and CH 4 
(CD 4 ) is 1000 times shorter than the time between collisions of 
two C 2 H radicals. The transverse flow arrangement in the cell 
allows high laser repetition rates with minimal photolysis of 
the same gas volume. With typical linear flow rates of 1.9 x 
10 20 molecules s -5 , about 90% of which is helium, the photolysis 
volume is replenished every 10 laser pulses. Kinetic experi- 
ments are performed at various methane (methane-*/) densities 
((0.2-2.5) x 10 17 cm -3 ). The total helium density is in the 
range (0.4-2.4) x 10 18 cm -3 , and the methane density is held 
at 0.7 x 10 17 cm“ 3 to test for a pressure dependence. In the 
experiments involving just the reaction C 2 H 4- C 2 H 2 , the 
acetylene number density was in the range (0.2-6.0) x 10 !5 
cm** 3 , and the helium number density was in the range (0.32 — 
3.2) x 10 ,8 cm~ 3 . For the highest acetylene pressures, the C 2 H 
concentration is estimated to be no greater than 6.7 x 10 ! 1 cm -3 . 
Therefore, a 300-9000-fold excess of C 2 H 2 with respect to C 2 H 
is always present. Contributions from radical— radical reactions 
can be neglected since the time between collisions for two C 2 H 
radicals is 1000 times longer at the highest C 2 H density than 
the time between collisions for C 2 H and C 2 H 2 . 

A high-resolution color center laser tuned to the Qn(9) line 
at 3593.68 cm -1 of the A 2 n— X 2 X transition probes the 
transient concentration of ethynyl radical in absorption. 9 A 
scanning Fabry-Perot spectrum analyzer is used to ensure that 
the color center is running on one longitudinal mode, and a 
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home-built scanning Michelson interferometer wavemeter 10 is 
used to monitor the color center’s wavelength. 

The probe beam, after three to five multipasses, is directed 
onto a 50-MHz 77 K GetAu detector which has a 20-mm 2 
sensitive area. The transient signals are amplified and then 
coadded using a 100-MHz digital oscilloscope. Typical single- 
shot traces have a signal-to-noise ratio of 20-30. For a typical 
run, transient signals from 1000 excimer pulses are averaged. 
The amplitude of the transient C 2 H signal is found to be linear 
with color center probe power. 

In all experiments, acetylene, helium, and methane or 
methane-d 4 were flowed through a mixing cell before entering 
the reaction cell. Helium is used to thermally equilibrate the 
mixture with the cell walls. In earlier experiments, sulfur 
hexafluoride, SF6, was used to vibrationally and electronically 
quench the C 2 H ; 3 6 this was omitted from this study because it 
was found to have no effect on the relatively slow ground-state 
removal rates, as discussed below. All gases are obtained 
commercially with the following purities: He, 99.99%; C 2 H 2 , 
99.6%; CILu 99.99%; CD 4 , 99.6%. The acetone in the C 2 H 2 is 
removed by passing the gas through an activated charcoal filter. 
Partial pressures of each gas are determined by calibrated mass 
flow meters and the measured total pressure inside the cell. With 
the use of isopentane as a cooling solvent, temperatures as low 
as 143 K can be reached. For measurements taken above 300 
K, heated water was used as the solvent. 

Analysis of Kinetic Data 

In this study we observe reactions from the ground state of 
C 2 H (X 2 Z + (0,0,0)) directly. For accurate measurements it is 
necessary that the electronic and vibrationally excited states of 
C 2 H be fully quenched. Work done by Glass et al. under similar 
conditions reports that relaxation of the C 2 H (X(0,0,1)) state 
occurs in approximately 1 /*s. 9 Therefore, if sufficient time has 
elapsed, complete vibrational and electronic relaxation should 
have occurred before any ground-state measurements are made. 
The data is fit only beginning with times a factor of 3 longer 
than the rise time to ensure complete relaxation of upper 
vibrational states of C 2 H. In previous studies SF6 was used as 
a vibrational quencher. 3 * 6 Temperature dependent measurements 
with and without SF6 were taken, and the measured rate 
constants were found to be equal. Therefore, SF 6 was not used 
for this study. 

Prior to measuring rate coefficients for the reaction C 2 H 4- 
CH 4 (CD 4 ), an accurate temperature dependence investigation 
of the rate constant for 

C 2 H + C 2 H 2 — C 4 H 2 + H (5) 

had to be completed. Earlier measurements could only be made 
down to 170 K, and here it was desired to extend the results to 
at least 143 K. 3 The experiments were done under pseudo- 
first-order conditions where [C 2 H 2 ] » [C 2 H] by a factor of 500- 
1500. The rate equation for reaction 5 integrates to 

1 [C : H], = [C 2 H] 0 exp(— * obs r) (6) 

where 

*obs = WC 2 H:] (7) 

The observed rate coefficients, /r 0 bs, are calculated by fitting 
the observed decay traces to a single-exponential decay plus a 
constant, eq 8, to fit the zero level of the base line 

>• = A exp (-k obs t) + constant 


( 8 ) 
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Figure 3. Plot of observed C : H removal coefficient *„ b s versus C : H : 
concentration. These data were taken at 155 K, and k c , H , for this plot 
is (1.6 ± 0.2) x lO' 10 cm 3 molecule ' s '. 


where A is the pre-exponential factor, t is the time, and k ob . is 
the observed rate coefficient; see Figure 2. The fits for k 0 bs are 
then plotted against their respective [C 2 H 2 ] concentrations, and 
the gradient is fco h 4 see Figure 3. The uncertainty in k Cz u : is 
calculated on the basis of the uncertainties in fitting k 0 bs , 
measuring [C 2 H 2 ], and measuring the temperature. This leads 
to relative errors in fcoHj of 10—13%, depending on the 
temperature. 

Table 1 is a summary of the rate constants measured ror L 2 H 
+ C 2 H 2 down to 143 K. The data from ref 3 and the new data 
can be fit approximately to the Arrhenius equation. However, 
the data is best fit by the equation kc : h : = 8.6 x 10 T 
exp[(474 ± 90 )/T] cm 3 molecule -1 s -1 ; see Figure 4. 

All experiments involving methane were performed under 
pseudo-first-order conditions in which [CRj], [CD 4 ], and [C 2 H 2 ] 
» [C 2 H]. The rate of change of [C 2 H] can be expressed as 

d[C 2 H]/dr = — [C 2 H](/r CXj [CX 4 ] + ^ c ,hJ C 2 H 23) (9) 

where CX 4 = CH 4 or CD 4 . After integration 

[C 2 H], = [C 2 H] 0 exp(-fc obs r) (60 


k obs = k c 4CX,) + k c ^[C 2 U 2 ] (10) 

Kbs ~ *C,H,l C 2 H 2] = ^CX 4 l CX 4l = k 'cX 4 ( 1 1 ) 

The observed decay rates, k 0 bs , are obtained by fitting the , 
observed traces to eq 8. They are then corrected for the 
contribution of the C 2 H 2 precursor reacting with C 2 H using both 
previous measurements and new measurements down to 143 
K, as noted in eq 1 1. Since C 2 H reacts 100 times faster with 
C 2 H 2 then it does with CX 4 , the contribution to k 0 bs from C 2 H 2 
4- C 2 H is kept below 40%. The values of k' cx 4 are plotted 
against their respective CX 4 concentrations, and a linear least 
squares fit is used to determine fccxji see Figure 5. The 
uncertainty in itcx 4 is calculated by combining the accumulated 
uncertainties in the corrected decay fits, which include errors 
in fcc, H „ with the uncertainties in temperature and the measure- 
ment of [CX 4 ]. The error in kcx 4 is typically 10-20%. 


TABLE 1: Summary of Rate Constants for C 2 H + C 2 H : 
from 143 to 359 K __ 

temp k c , H, temp , 

(K) (cm 3 molecule -1 s -1 ) (K) (cm 3 molecule ' 1 s ) 

359" (1.3 ±0.1) x 10" 10 

350* (1.2 ± 0.2) x 10 -l ° 

296" ( 1 .4 ± 0. 1 ) x 10 -10 

295* (1.3 4 0.2) x 10 -10 

273* (1.2 ±0.2) x 10 -10 

263* (1.2 4 0.2) x 10" 10 

256* (1.3 4 0.2) x 10-'° 

245* (1.3 4 0.2) x 10 -10 

242* (1.2 4 0.2) x 10-'° 

239* (1.3 4 0.2) x 10"'° 

233* (1.2 4 0.2) x 10 -10 

232* (1.3 4 0.2) x lO" 10 

225* (1.2 4 0.2) x 10"'° 

223* (1.3 4 0.2) x 10-'° 

223* (1.2 4 0.2) x 10 10 

213* (1.4 4 0.2) x 10 -10 

« Measurements in this work. * Measurements made in a previous 
study . 3 



Table 2 summarizes the kinetic measurements in these 
experiments for ken, and kcui from 154 to 359 K. The date 
indicate that the rate coefficients ken, and k C o. have a positive 
temperature dependence from 154 to 359 K, which can be 
expressed as Jfccu, = 0-2 ± 0.1) x 10 - " exp[(-491 ± 12)/T 
and ^cd 4 = (8.7 ± 1.8) x 10 -12 exp[(-650 4 61)/21 cm 
molecule -1 s -1 , respectively; see Figure 6. The reaction of C 2 F 
+ CD 4 exhibits a significant kinetic isotope effect at 300 K o 
kcnJkcDi = 2.5 ± 0.2. 

Discussion 

Previous measurements of the reaction C 2 H + CHa wen 
performed at room temperature and are summarized in Tab! 
3. Rendlund et al. monitored the CH (A-X) product chemilu 
minescence from the C 2 H 4 0 2 reaction with and without CH4. 
They measured a room temperature value of kcn 4 = (4.8 ± 1.0 
x 10" 12 cm 3 molecule -1 s -1 . Okabe photolyzed C 2 H 2 at 14 
nm and determined a ratio of k C Hjkc : H : = 0.032 ± 0.001 J 
Using the present value of kc z h ; at room temperature of ( 1 .3 : 
0 2) x 10“ 10 cm 3 molecule -1 s’ -1 , one obtains a value for key 
of 4.2 x 10 -12 cm 3 molecule -1 s -1 . 12 Laufer photolyz 


2ir 

(1.3 ±0.2) x 1U 

207* 

(1.3 4 0.2) x 10 -l ° 

201* 

(1.3 4 0.2) x 10 -10 

199* 

(1.3 4 0.2) x 10 -10 

197* 

(1.4 4 0.1) x 10 -10 

193* 

(1.3 4 0.2) x lO -10 

193* 

(1.3 4 0.2) x 10 -10 

191“ 

(1.5 4 0.1) x 10' 10 

188* 

(1.2-4 0.2) x 10 -10 

185* 

(1.4 4 0.2) x 10 -10 

170* 

(1.4 4 0.2) x 10 -l ° 

163" 

(1.8 4 0.2) x 10 -10 

155" 

(1.6 4 0.2) x 10 -10 

153* 

(1.8 4 0.2) x 10 -10 

143" 

(1.9 4 0.2) x 10 -10 
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61)/7] cm 3 molecule' 1 s' 1 , respectively. 


TABLE 2: Summary of Rate Constants for C 2 H + CHt 

(CD 4 ) from 154 to 359 K 


temp 

(K) 

*ch 4 

(cm- molecule" 1 s" 1 ) 

temp 

(K) 

*CD; 

(cm 3 molecule" 1 s" 1 ) 

359 

(3.1 ±0.4) x 10' 12 

359 

(1.5 ±0.2) x 10' 12 

300 

(2.3 ±0.1) x 10' 12 

300 

(0.9 ±0.1) x 10' 12 

226 

(1.2 ±0.1) x 10' 12 

227 

(5.2 ±0.1) x 10' 13 

190 

(8.7 ± 0.2) x 10' 13 

190 

(2.8 ± 0.6) x 10' 13 

159 

(5.2 ± 0.6) x 10' 13 



157 

(4.7 ± 1.0) x 10' 13 



154 

(5.5 ± 1.0) x 10" 13 



CF 3 C 2 H in the presence of CH 4 , 

and the 

product [C 2 H 2 ] was 


monitored by absorption at 152 nm over time . 13 Laufer obtained 
a item = (1.2 ± 0.2) x 10 ' 12 cm 3 molecule ' 1 s' 1 . Finally, 
Lander et al. used a diode laser to measure the transient 
depletion of C 2 H (X 2 Z+ (0,0,0)) in absorption . 14 They measured 
a fern, = (3.0 ± 0.3) x 10' 1: cm 3 molecule ' 1 s " 1 at 300 K. 

The present rate coefficient measurements for CjH + CH 4 
at 300 K and from 154 to 359 K are &ch 4 = (2-3 ± 0.1) x 
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TABLE 3: Summary of Rate Constants for C 2 H + CH 4 at 

Room Temperature 


ref 

kctii (cm- molecule" 1 s" 1 ) 

temp (K) 

this work 

(2.3 ±0.1) x 10' 12 

300 

Renlund 11 

(4.8 ± 1.0) x 10' 12 

298 

Laufer 13 

(1.2 ±0.2) x 10' 12 

297 

Okabe 12 

4.2 ± 10' 12a 

298 

Lander 14 

(3.0 ±0.3) x 10' 12 

298 


" Recalculated using the measured ratio and the present-day C:H + 
C 2 H 2 rate constant (1.3 ± 0.2) x 10" 10 cm 3 molecule' 1 s' 1 . 


10' 12 and ken* = (1-2 ± 0 . 1 ) x 10' 11 exp[(-491 ± 12)/7J 
cm 3 molecule ' 1 s' 1 , respectively. From the Arrhenius fit, the 
energy of activation for C 2 H + CH,; is equal to 4.3 ±0.1 kJ 
mol' 1 . Rate coefficient measurements for C 2 H + CD 4 at 300 
K and from 227 to 359 K result in a k c d 4 = (0.9 ± 0.1) x 
10 ' 12 and kcDi = (8.7 ± 1.8) x 10 ' 12 exp[(-650 ± 61)7] cm 3 
molecule ' 1 s' 1 , respectively. The energy of activation for C 2 H 
+ CD 4 from 227 to 359 K is equal to 5.4 ± 0.5 kJ mol' 1 . No 
pressure dependence was found for &ch 4 at 300 K over the 
experimental range (0.4-2.4) x 10 18 cm ' 3 variation in helium 
buffer gas (15-75 Torr) and for a methane density of 0.7 x 
10 17 cm" 3 . Over this number density k C H 4 = (2-4 ± 0.2) x 
10 " 12 cm 3 molecules ' 1 s' 1 . Both of these reactions, C 2 H + 
CHa (CD 4 ), show a positive temperature dependence, consistent 
with reactions involving a hydrogen abstraction as the rate- 
determining step and a positive activation energy. 

Our results are most consistent with those of Lander et al. 
Both experiments monitor the direct disappearance of C 2 H. The 
difference in k c h 4 at 300 K between this work and the work by 
Lander et al. can be attributed to experimental error in 
determining number density and error in fitting the observed 
signal. In the experiment by Renlund et al. the chemilumines- 
cence of CH (A 2 A) produced by the reaction C 2 H + 0 2 was 
used to monitor [C 2 H]. They were able to measure k CH , by 
keeping the concentration of 0 2 constant while varying the 
concentration of CHa. Renlund et al. measured a value for k c h 4 
at 300 K that is a factor of 2 greater than the rate constant 
reported in this work. They reported in a later publication that 
the reason for the larger £chj was due to vibrationally and or 
electronically excited C 2 H producing the CH (A 2 A-X 2 I1) 
emission . 15 Okabe measured /cchAc ; H; to be 0.032 ± 0.0018. 
This leads to a value of fc C H 4 that is ~2 times faster at room 
temperature than the ken, reported in this work when using (1.3 
± 0.2) x 10 ~ 10 cm 3 molecule ' 1 s ' 1 for kc : n v However, 
Okabe’s result involved measuring d>°c 4 H ; /0c 4 H : , which is the 
quantum yield ratio of diacetylene without (d>°c 4 H : ) and with 
(0 c 4 h 0 CRt in the mixture and not a direct measurement of the 
reaction of C 2 H with CH 4 . 

The work presented here on C 2 H + C 2 H 2 further supports 
the reasoning that C. 1 H 3 * is formed via an addition mechanism 
with no barrier to formation. The short-lived intermediate then 
immediately dissociates to C 4 H 2 4- H. A previous study on 
C 2 H + C 2 H 2 by this lab reported that within experimental error 
there was no temperature dependence over 170-350 K . 3 More 
recent results by Van Look et al. also showed no evidence for 
a temperature dependence from 295 to 450 K . 16 They measured 
rate coefficients of the reaction C 2 H + 0 2 by monitoring CH- 
(A 2 A — X 2 n) chemiluminescence. The rate coefficients for 
C 2 H with C 2 H 2 were taken from the y-intercepts of their k' 
versus [0 2 ] plots. An Arrhenius fit to their data is reported as 
fc c , H , = (1.3 ± 0.2) x 10 ' 10 exp[(0 ± 10)/7] cm 3 molecule ' 1 
s' 1 . Despite previous measurements for the rate coefficients 
of C 2 H + C 2 H 2 , which show no definite temperature dependence 
above room temperature , 3 - 916 - 17 new data presented in this study 
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does reveal a small negative temperature dependence in the rate 
coefficients for CSH 4- C 2 H 2 . 

The measured kinetic isotope effect at 300 K is IcchAcDi = 
2.5 ± 0.2. Unfortunately, no transition-state calculations have 
been done on this reaction to be able to compare theory with 
experiment. Judging by the magnitude of the isotope effect and 
the reaction mechanism involved, the reaction C 2 H -4- CH 4 (CD 4 ) 
exhibits a large primary isotope effect in which the main 
contribution is the difference in zero-point energies between 
the initial states and the transition states . 18 " 20 Hopefully this 
research will encourage theoretical studies of C 2 H + CH4 (CD 4 ) 
for comparison. 

In addition to the kinetic isotope effect, there is the possibility 
for tunneling to occur in this reaction between the H (D) atom 
of methane and the carbon atom of the C 2 H radical. Unfortu- 
nately, no curvature is detected in the Arrhenius plot over 1 54 
359 K. This does not eliminate the possibility of tunneling 
taking place. If the experiment could be done over a larger 
temperature range, the possibility of detecting curvature in the 
Arrhenius plot would be more favorable. 

Results based on the effect of reaction 2 on photochemical 
models of Titan are forthcoming. There is little question that 
this reaction plays a key role in determining the ethane 
concentration on Titan, but recent experiments by Mordaunt et 
al? ] have uncovered a new CH4 photolysis channel that will 
affect current photochemical models. Mordaunt et al. have 
shown that the direct photodissociation of CH4 to CH3 + H is 
the main source of methyl radicals in Titan’s atmosphere. Both 
the discovery by Mordaunt et al. and the data in this work should 
help to clarify the production of ethane in Titan’s atmosphere. 
Updated versions 2 of the Yung et al 1 model and more low- 
temperature rate coefficients of pertinent reactions are needed 
to determine whether a complete analysis will match forthcom- 
ing spacecraft observations. 
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